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Abstract. The final state of Au + Au collisions at -^s = 130 AGeV at 
RHIC has been reconstructed within the framework of the Buda-Lund hydro- 
dynamical model, by performing a simultaneous fit to final data on two-particle 
Bose-Einstein correlations of the STAR and PHENIX Collaborations, and fi- 
nal identified single particle spectra as measured by the PHENIX Collabo- 
ration. The results indicate a strongly three dimensional expansion, with a 
four-velocity field that is almost a spherically symmetric Hubble flow. We find 
large transverse geometrical source sizes, Rq = 9.8 ±1.2 fm, relatively short 
mean freeze-out time, tq — 6.1 ± 0.3 fm/c and a short duration of particle 
emission, Ar — 0.02 ± 1.5 fm/c. Most strikingly, we flnd an indication for 
a hot central part of the hydrodynamically evolving core, characterized by a 
central temperature of To = 202 ± 13 MeV, that is close to (or even above) the 
deconfinement temperature of the quark-hadron phase transition. The best fit 
indicates a cold surface temperature of Tj = 110 ± 16 MeV. When the pos- 
sibility of the hot center is excluded, the confidence level of the fit decreases 
from 28.9% to 1.0 %. Predictions are made for the rapidity dependence of the 
slope parameters and for the transverse mass dependence of the rapidity width 
of the single particle spectra, and the transverse velocity dependence of the 
non-identical particle correlations. 
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1. Introduction 

The reconstruction of hadronic final state from the measured single particle spectra 
and two-particle correlation functions is of great current research interest in high 
energy heavy ion collisions. The goal of these studies is to characterize the final state 
with the help of a few simple parameters like the mean freeze-out temperature or the 
mean transverse flow, or the geometrical size of the system at the mean freeze-out 
time. From the properties of the hadronic final states, one aims to identify one (or 
more) new phases of hot and dense hadronic matter in the collisions of the biggest 
nuclei at the highest available bombarding energies. 

Space-time picture reconstruction of particle emitting sources in high energy 
physics is based on improved methods of intensity interferometry, a technique in- 
vented originally by the radio astronomers R. Hanbury Brown and R. Q. Twiss [ 0| 
to measure the angular diameter of main sequence stars. Intensity correlations in 
high energy physics are measured in momentum space, in contrast to stellar inten- 
sity interferometry, where the correlations are determined in the coordinate space. 
Intensity correlations of identical particles appear due to quantum statistics as well 
as Coulomb and strong final state interactions. Dominated by quantum statistical 
effects, correlations of identical bosons are frequently referred to as Bose-Einstein 
correlations (and the name of fermionic correlations is Fermi-Dirac correlations). 
The radius parameters of the two-particle Bose-Einstein correlation functions are 
frequently referred to as HBT radii to honor Hanbury Brown and Twiss. In particle 
physics, correlations of pions with small opening angles were observed first by G. 
Goldhaber, S. Goldhaber, A. Pais and W. Lee, [ || hence they are also referred to 
as GGLP correlations. 

Our tool for the space-time picture reconstruction is the Buda-Lund hydro (BL- 
H) model, introduced in refs. [ ^, |^ and reviewed recently in ref. [ ^. See also refs. [ 
0, |, I for recent reviews on particle interferometry in high energy physics. 

First of all, let us remind the readers about the difficulty of the interpretation 
of the experimental data on correlations and spectra in high energy heavy ion colli- 
sions. High expectations seem to exist each time, when a new accelerator starts its 
data taking. Our hopes suggest that we enter a whole new land to explore, and we 
tend to forget about the landscape that we have just left behind. 

Let us remember, that the "RHIC HBT puzzle" [||] is in fact not a RHIC specific 
phenomena: a similar discrepancy between predictions and final data prevailed 
already at CERN SPS, see refs. Essentially, the RHIC HBT puzzle is that 

the ratio of the outward and the sideward HBT radii was measured to be Ro/Rs ~ 
1.0, see refs. [ |l^, and this result was in contrast to certain theoretical 
expectations. The "CERN SPS HBT puzzle" could be similarly formulated: Why 
the difference between the outward and the sideward HBT radii is zero at CERN 
SPS? This question can be based on Fig. 4 of ref. [ see ref. [ 111] for the most 
recent data with such a behavior in Pb-|-Pb collisions at CERN SPS, valid not only 
for pions but also for kaons. At CERN SPS, the situation seems to be strikingly 
similar to the happenings at RHIC. Perhaps the story at CERN SPS is so old by 
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now, that its moral has been almost forgotten. 

Apparently, failed predictions of the RHIC HBT data can be assigned to models 
that were not tuned to successfully describe the single particle spectra and the 
two-particle Bose-Einstein correlation functions at CERN SPS. On the other hand, 
models that worked well at CERN SPS describe RHIC spectra and HBT data rather 
well, see e.g. refs. [ |3 |l| . 

There were two classes of predictions for the measurable HBT radius parameters 
in Au + Au collisions at RHIC. The predominant expectation was, that a soft, 
long lived, evaporative quark-gluon plasma phase will be produced, and the large 
duration of the particle emission can be observed from the big increase of the out 
component of the Bose-Einstein correlation functions, i?out S> i?sidci a- phenomena 
first observed by S. Pratt in refs. [ Similarly motivated calculations, with 

more realistic geometry and initial conditions were performed by Bertsch [ ^ and 
collaborators as well as by Gyulassy and Rischke [ ^ . 

However, there was another, less well known class of predictions for RHIC: 
instead of predicting i?out/-Rsidc ^ 1, refs. [ ||, ^ predicted i?out/-Rsidc ~ 1 ■ 
It is also mentioned there how the sudden freeze-out of hadrons is related to an 
explosive particle production from a supercooled quark-gluon plasma [ |2^, or a 
quark matter [ 27|, where the gluonic degrees of freedom are not active. Such a 
picture may emerge from a quasi-particle picture of a QGP, where the quarks and 
the gluons dress up in the vicinity of the phase transition, the constituent quark 
mass being of the order of 300 MeV, with a gluon mass nig ^ T^, [ The 
sudden, explosive particle production, the hard equation of state and the quark 
combinatorics of particle yields [ ^ supports such a scenario at RHIC. 

As the characteristic nucleation times of hadronic bubbles inside a quark gluon 
plasma are of the order of 100 fm/c, which is an order of magnitude larger than 
the characteristic life-time of the expanding system, bubble formation is not fast 
enough to keep the system close to the Maxwell construction and near-equilibrium 
phase transition. Instead, a negative pressure state developes very soon which then 
decays due to its mechanical instability, the cavitation. Such process may happen 
through a time-like deflagration and a process was predicted to end in a pion flash, 
with a short, 1-3 fm/c duration of particle emission. See ref. [ ^ for greater details 
and signatures of this process. 

In this manuscript, we present an indication for the existence of a transient 
quark matter state at RHICl. 

The structure of the body of the paper is as follows. In section ||, we highlight 
the most important fitting formulas of the BL- hydro model. In section H, the fits to 
the single particle spectra and the two-particle Bose-Einstein correlation fimctions 
are shown. In section ^ we summarize the results, and discuss their interpretation of 
the results and the correlations between the various fitting parameters. In section ^ 
we predict the rapidity dependence of the observables. We also make a prediction for 
the transverse mass, and rapidity dependence of non-identical particle correlation 
functions. Finally, we conclude. 
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2. Single particle spectra and two particle correlations from 
the Buda-Lund hydro model 

Our direct aim is to reconstruct the hadronic final state from the measurable single- 
particle spectra and two-particle correlation functions. From this reconstructed 
final state and the knowledge of the equation of state of hot and dense hadronic 
matter (e.g. from lattice QCD calculations) one can, in principle, reconstruct the 
initial state of the reaction by running the (relativistic) hydrodynamical equations 
backwards in time, and determine if this initial state had been in the QGP phase or 
not. Here we report on such a direct reconstruction of the hadronic final state within 
the framework of the Buda-Lund hydro (BL-H) model, but we do not consider the 
more indirect reconstruction of the initial state. 

The BL-H model is a hydrodynamical parameterization of the hadronic final 
state, that has to be clearly distinguished from a fully developed, time dependent so- 
lution of relativistic hydrodynamics. However, the BL inspired a whole new series of 
non-relativistic and relativistic, simple analytic solutions of fireball hydrodynamics. 
The non-relativistic hydro solutions correspond to the non-relativistic limit of the 
BL hydro model, with time-dependent model parameters, ^ |3^, 

It is very interesting to note, that the governing equations for the scale parameters 
are similar to that of the Zimanyi-Bondorf-Garpman hydrodynamical solution [ |3^ 
and its ellipsoidally symmetric generalization [ |3^ . 

However, in the relativistic domain, only the coasting (accelerationless) hydro 
solutions were found analytically until now [ |3^, ^ , in an attempt to figure 
out the governing equations for the BL-H model parameters. Although Bjorken's 
well known solution of relativistic hydrodynamics [ ^ belongs also to this class of 
accelerationless solutions, and the new coasting relativistic solutions include finite 
1-fl dimensional solutions [Js^, 1-1-3 dimensional solutions with cylindrical [ |39|| 
and ellipsoidal symmetry [ [40| , the search is still going on for even more realistic, 
relativistic hydro solutions that can interpolate from the non-relativistic domain to 
the relativistic one even if the acceleration of the matter is significant. 

Furthermore, the Buda-Lund flow profile, with a time-dependent radius param- 
eter Rq, was recently shown to be an exact solution of relativistic hydrodynamics 



of a perfect fluid at a vanishing speed of sound [ 42 . It turned out [ ^ ^ , that 
the flow field is a generalized Hubble flow and the average transverse flow at the 
geometrical radius is formally similar to Hubble's constant that characterizes the 
rate of expansion in our Universe, (ut) ~ ^tRc — ItH [ This emphasizes the 
similarity between the Big Bang of our Universe and the Little Bangs of heavy ion 
collisions. 

The invariant single particle spectrum is obtained [ ^ ^ from BL-H as 

(^fl g t 

^'^^^ = 27rmtdmt dy ^ (2^ ^ ^ ^ ^^^^ f ^^(S)fc^ ' 



T{x) T{x 

where all the terms have an intuitive, but mathematically well defined meaning 
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see ref. [ |^ for more detailed definitions. The invariant form of the two-particle 
Bose-Einstein correlation function (HBT) of the BL-H is found in the binary source 
formalism [ ^ |^ as: 

C2(ki,k2) = 1 + A,0(g||) exp(-Qfi?J-Qii?L-Qii?l) . (2) 

The dependence of the fit parameters on the value of the mean momentum of the 
pair is suppressed, see ref. [ ^ for the complete set of definitions of the variables 
and the radius parameters in terms of the BL-H model parameters. The pre-factor 
ri(Q||) of the BECF induces oscillations within the Gaussian envelope as a function 
of Q\\. This oscillating pre-factor satisfies < ^{Q\\) < 1 and ri(0) = 1. In 
practice, the period of oscillations is larger, than the corresponding Gaussian radius, 
so the oscillations are difficult to resolve. The above BL-H form of the two-particle 
correlation is explicitly boost-invariant, as all the relative momentum dependent 
variables Q\\, Q=, Q± and all the corresponding radius parameters, R=, R± 
are defined in an explicitely boost invariant manner. The BL correlation function 
can be equivalently expressed in the frequently used, but not invariant Bertsch-Pratt 
(BP) form in the LCMS frame [ within the 17 = 1 approximation: 

C2(ki,k2) = 1 + A, exp [-i?2Q2 _ i?2g2 _ ^2q2 _ 2i?2jQ^Ql] . (3) 

The above formulas for the BECF and IMD, as were used in the fits, have been 
introduced in refs. ^, |4^, and summarized recently in ref. [ ||. The analytic 
formulas that relate the BL-H model parameters to the above forms for the spectra 
and correlation functions, are given by eqs. (84-105), (115-118) and (129-140) of 
ref. [|. Note, however, that eq. (132) of ref . [|| contains an unfortunate misprint, 
{{ut) + (AT/T),.) in the denominator should be replaced by ((ut)^ + (AT/r)^), so 
the correct form of eq. (132) of ref. [ |^ given explicitely by eq. (||) of the present 
manuscript. 



3. Buda-Lund fits to RHIC-1 spectra and correlations 

Here, we reconstruct the space-time picture of particle emission in Au + Au colli- 
sions at RHIC within the BL-H framework, by fitting simultaneously the PHENIX 
and STAR final data on two-particle correlations and single-particle spectra pre- 
sented in refs. [ ^ |lj, |l2|. The BL model, in certain domain of the parameter 
space [ ^, features a scaling limiting behavior of all the HBT radius parameters, 

i?out — -Rsidc — -Riong ~ tq^^^. A uuiquc feature of the BL model is that it has 
been successfully tested against a detailed description of the single-particle spectra 
and the two-pion Bose-Einstein correlation functions in both h + p and Pb + Pb 
collisions at CERN SPS bombarding energies [ [l^ . 

The measured and the calculated single particle spectra are connected with the 
help of a core-halo correction factor cx l/\/X^ , where the experimental values of the 
intercept parameter A*(?/, mt) are to be taken from the measurements. In the lack 
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of these A*(y, mt) values we have utihzed their average A* for a core-hafo correctfon 
when fitting the PHENIX K~ and p~ spectra. In particular, the following average 
values were used for the various particle types: \^,{K) — 0.80 (estimated from 
the NA44 data on kaon-kaon correlations at CERN SPS [ 0, A*(p) = 0.995 (the 
fraction of long lived resonances that decay to anti-protons is neglected). For pious, 
we have utilized the A*(mt) values given in refs. [ 



3.1. Improving on earlier results 

Note also that we have performed the data analysis within a Gaussian approxima- 
tion to the oscillating prefactor, improving on our earlier results [ ^ , where we 
have utilized the 17 = 1 approximation. As the oscillating prefactor Q depends only 
on (3||, a Gaussian approximation to Q results in a a Gaussian form of the Bose- 
Einstein correlation function in the Buda-Lund variables, but with a modification 
of the invariant longitudinal size of the source, 



1 



A77 



cosh (rj) 



(4) 



where is given by eq. (138) of ref. This modified invariant longitudinal radius 
2 2 

parameter, q replaces in eqs. (116-118), when making the transformation 
from the Buda-Lund radius parameters to the experimentally determined Bertsch- 
Pratt radius parameters following the lines of ref. [ ^ . 

Furthermore, we improved on our earlier results [ ^ by taking into account an 
mt dependent core-halo correction for the PHENIX spectra and correlations, and by 
fitting the absolute normalization of the single particle spectra in both experiments, 
properly utilizing the degeneracy, fugacity and quantum statistical factors. This 
allows us to extract the chemical potential in the center of the fireball, in contrast 
to our earlier fits [ ^ where the absolute normalization of the particle spectra and 
the central value of the chemical potential distribution were not yet determined. 
In ref. [ we have attempted to determine the central value of the chemical 
potentials from the preliminary data. We improve also on this analysis by using the 
final, published data, by releasing the central value of the pion chemical potential, 
that was previously fixed to 0, and by using the correct value of the g ^ 2s + I the 
degeneracy factors, i.e. g — (1, 1,2) for (tt", K~, p~). 

Fig. 1 illustrates the best combined fit to the single particle spectra of negative 
pions, kaons and anti-protons of PHENIX, as well as to the transverse mass depen- 
dent HBT radii of STAR and PHENIX. Thus the hypothesis, that pions, kaons and 
protons are emitted from the same hydrodynamical source is in a good agreement 
with the fitted data, and the PHENIX and the STAR datasets are compatible with 
one another. The parameters of the combined fit to PHENIX and STAR data are 
summarized in Table 1. 

A very important improvement over our earlier results is that now we have 
utilized the final PHENIX and STAR data points, in contrast to the earlier results, 
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BL hydro 


STAR+PHENIX 


STAR+PHENIX 


STAR+PHENIX 


T 

J.. 


TT 


TTT 


Rg [fm] 
To [&n/c] 
At [fm/c] 
A?7 

(~T^)l" 

/AT\ 
\ — /t 


1 08 ± 17 
9.8 ± 1.2 
6.1 ± 0.3 

0.02 ± 1.5 
2.5 fixed 

0.84 ± 0.24 
1.3 ± 0.4 


80 ± 12 
7.8 ± 0.8 
7.6 ± 0.2 

0.10 ± 0.9 
2.5 fixed 

0.10 ± 0.04 
0.4 ± 0.1 


1 1 n fivprl 

J. J. U llXtJl.! 

69 ±0 11 
7.3 ± 0.8 
8.9 ± 0.2 

0.07 ± 0.6 
2.5 fixed 

0.01 ± 0.01 
0.1 ± 0.1 


Mr [MeV] 
(4' [MoV] 

M?; [-McV] 


75 ± 19 
107 ± 14 

305 ± 41 


104 ± 14 
172 ± 17 

412 ± 17 


118 ± 12 
206 + 12 

437 ± 17 


xVNDF 


74/68 = 1.09 


99/69 = 1.43 


127/69 = 1.84 


CL 


28.9% 


1.0% 


2.7e-3% 



Table 1. Source parameters from simultaneous fits of final Au + Au RHIC data 
of PHENIX and STAR on particle spectra and HBT radius parameters with the 
Buda-Lund hydro dynamical model. 

where the preliminary data points wer(^ utilized. Unique minima are found and a 
statistically acceptable x^/NDF is obtained for the combined fit of both PHENIX 
and STAR datasets. Furthermore, the /NDF values are good not only on the 
combined data set, but on each of the fitted spectra and HBT radii, perhaps with 
the exception of the out radius parameter R^ at the largest tranverse mass value, 
where the fit overestimates the measured point by about 2.5 standard deviations. 

The improvements did lead to a refinement of the fitted paranic;ters and their 
values, but, typically, the changes were within 3 standard deviations of the errors, 
with the exception of the central value of the freeze-out temperature, which in- 
creased significantly, from the preliminary RHIC average of Tq 142 ± 4 MeV to 
the final Tq = 202 ± 13 MeV, which is a significant, close to 5 standard deviation 
modification. Due to the above improvements, the analysis of the final STAR and 
PHENIX data now indicates the existence of a hot central part of the core, which 
was not resolved in the preliminary analysis. Clearly, more data points would be 
very useful to confirm or invalidate this observation. In particular, HBT radius pa- 
rameters of pions and kaons in the « 1 GeV domain would provide very useful 
constraints on the parameters of the BL-hydro model. 

4. Discussion 

We emphasize that the way how the data were presented are not ideal for fitting the 
BL-hydro model parameters. In particular, the BL-hydro model parameters can be 
determined easier, if the rapidity dependent transverse mass spectra are given for as 
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BL hydro 


Best fit 


To 140 MeV 


To = 110 MeV 


parameters 


Surface values 


Surface values 


Surface values 


Ts [MeV] 


110 ± 16 


127 ± 5 


108 ± 1 




0.73 ± 0.06 


0.62 ± 0.04 


0.57 ± 0.03 


[MeV] 


-25 ± 21 


34 ± 15 


63 ± 12 


[MeV] 


6 ±29 


103 ± 16 


152 ± 12 


M? [MeV] 


204 ± 49 


343 ± 17 


402 ± 17 


Conf. level 


28.9 % 


1.0 % 


2.7 10-^% 



Table 2. Calculated parameters, the surface values of the temperature, flow and 
chemical potential distributions, as evaluated from the simultaneous fits of Table 
1 to final Au + Au RHIC data of PHENIX and STAR on particle spectra and 
HBT radius parameters with the Buda-Lund hydrodynamical model. The surface 
temperature is Ts = 70/(1 + (AT/T)^), the average surface three-velocity is (/3t) = 
{ut) I a/I -I- and the surface chemical potentials are given by \i,s = I^q — To/2 for 
all particles. 

many particles, as possible, as the BL model contains terms that explicitely break 
the longitudinal boost invariance of the source. For an example of such an analysis, 
see ref. [ |l^ . The BL HBT radii also contain small correction terms that depend on 
the deviation from mid-rapidity in an explicite manner. The violation of complete 
longitudinal boost invariance with a finite space-time rapidity distribution is in fact 
an important feature of the BL-hydro model, that can be utilized to distinquish the 
BL-hydro from longitudinally boost invariant sources, like Bjoken's hydro solution, 
or, the blast- wave model of ref. [Q. Note that the terminology "blast-wave" was 
introduced in ref. [ ^ , which generalized phenomcnologically the Zimanyi-Bondorf- 
Garpman (ZBG) solution [ ^ of non-relativistic hydrodynamics to relativistic flow 
profiles. It is very interesting to observe this as the BL-hydro model also corresponds 
to a relativistic extension of the ZBG solution. 

The major difference between the final state of heavy ion collisions at RHIC and 
at CERN SPS seems to be not only the increased freeze-out time and the increased 
transverse fiow or Hubble constant and relevantly larger tranverse geometrical ra- 
dius at RHIC, but, most strikingly, the existence of a hot center located close to 
the beam axis, which evaporates particles with a temperature that is very close to, 
or above the deconfinement temperature, Tq = 202 ± 13 MeV. This effect was not 
seen at CERN SPS, when the NA44, NA49 and WA98 data on particle spectra and 
correlations were analyzed with the help of the BL-hydro model [ |T^. We have 
checked, that this minimum satisties eqs. (16-18) of ref. [ 0, i.e. the conditions 
for the validity of the saddle-point approximation are satisfied. 

We tried to determine the significance of this result by setting the central tem- 
perature to the value of To = 140 MeV, artificially requiring that the temperature 
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of the hottest zone at RHIC be similar to the value of Tq, that were found in h + p 
and Pb + Pb reactions at CERN SPS [ |l|, A reasonable fit was obtained, how- 
ever, the confidence level of the fit decreased from 29 % to 1.2 %, and the /NDF 
increased correspondingly, see the second column of Table 2. The blast- wave model 
describes the preliminary transverse mass dependence of the STAR data on neg- 
ative pions, kaons, anti-protons and anti-lambdas with Tq = 120^25 MeV, and 
(A) = 0.52 

iao8 [ HI- Within errors, we recover this minimum from the combined 
fits to the PHENIX final spectra and the final STAR and PHENIX HBT radius 
parameters, if we require a vanishing transverse temperature gradient in the BL- 
hydro model, or, if we fix by hand the central value of the freeze-out temperature 
distribution to Tq = 110 MeV, see the third column of Table 1. However, we find 
that this minimum is statistically not acceptable, as the confidence level of the fit is 
very small. We have done another check, when determining the stability of the fit 
results for the variation of the central value of the freeze-out temperature. We have 
investigated the stability of the surface temperature for the variation of Tq (best 
fit value versus Tq = 140 MeV fixed versus Tq = 110 MeV fixed). Interestingly, we 
found that the transverse temperature gradient, (AT/T)^ is strongly correlated to 
To, so a higher central temperature in the fit implies a larger temperature gradi- 
ent, in such a way, that the surface temperature, Tg — 7o/(l + {AT/T)r) remains 
within errors the same. We have also calculated the 3-velocity of the matter on the 
surface in the transverse direction, {(3t) = {ut)/\/l + {u)'j for a comparison, and 
the chemical potentials on the vt — Rq surface, that are given by /i^ ^ fiQ — Tq/2 
for all particles. The freeze-out temperature and the average transverse velocity 
at rt = Rg are remarkably stable parameters of the fits, but, only the results in 
the first column of Table 2 are statistically acceptable. According to this column, 
the chemical potential of negative pions and kaons approximately vanishes at the 
rt — Rg surface, and the chemical potential for anti-protons is significantly bigger, 
than zero on the surface. 

Furthermore, the surface temperature and the surface three- velocity was found 
to be similar to the average freeze-out temperatures that were obtained from the 



Regensburg model in ref. [ 50 1 and the blast-wave model analysis of the RHIC final 



state in ref. [ pll. Our results are qualitatively as well as quantitatively similar to 



the findings of Florkowski and Broniowski [ 20 , who analyzed the PHENIX and 
STAR single particle spectra at RHIC in a hydro model that includes a spherically 
symmetric Hubble flow and similar in spirit to the BL-hydro and the blast-wave 
models. They find that a value of To = 165 ± 7 MeV, ^f, = 41 ± 5 MeV and an 
average transverse flow velocity of (/3t) = 0.49 describes the STAR and PHENIX 
preliminary data on the single particle spectra not only for negative pions, kaons 
and anti-protons, but also for mesons and A-s, A-s, and K*-s, [ 

We have found a non- vanishing chemical potential for negative pions, kaons 
and anti-protons in the center of the fireball from the absolutely normalized single- 
particle spectra. The pion and kaon data were well described in all cases with a 
chemical potential that (within errors) vanishes on the rt = Rg surface of the fire- 
ball. These values together with the inhomogeneous chemical potential distribution 
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of eq. (124) of ref. [ |^ indicate a clear deviation from chemical equilibrium in the 
hadronic final state, as reconstructed within the Buda-Lund hydro model. 

The similarities and the differences between an effective Quark Matter (QM) 
stage and a Quark Gluon Plasma (QGP) phase have been summarized recently in 
ref. [ ^ . The observed short duration of particle emission and the large transverse 
flow at RHIC contradicts to the picture of a soft, long-lived, evaporative Quark 
Gluon Plasma phase, that would consist of massless quarks and gluons. However, 
the final state does not exclude a transient, explosive, suddenly hadronizing Quark 
Matter phase, that could be characterized by massive valence quarks, the lack of 
gluons as effective degrees of freedom, and a hard equation of state. 



5. Buda-Lund Predictions - What Next? 

In the above fits, we have utilized only mid-rapidity data points as given by the 
PHENIX and STAR collaborations. Additional measurements provide more strin- 
gent restrictions on the value of the fit parameters. 



5.1. (Pseudo-)rapidity dependent measurements 

Recently, BRAHMS has published the pseudorapidity distribution of charged par- 
ticles [ |5^, which helps to restrict the value of the width parameter Arj, and also 
to get a more precise handle on the difference between the central and the surface 
temperature of the fireball, as the broadening of the pseudo-rapidity distribution 
depends on this parameter. Preliminary results indicate an agreement with the 
value of A77 = 2.5 utilized in the fits presented here. 

The Buda-Lund hydro model predicts a specific a coupling between the rapidity 
and the transverse mass dependent single particle spectra . In particular, due to 
the finite longitudinal size of the expanding fire-tube it predicts ^ |^, that the 
effective slope parameter decreases in the target and the projectile fragmentation 
regions as 

Tcsiy) = — — 7^, (5) 

where (transverse mass dependent) slope parameter at mid-rapidity is given by 

= To + nit (ut) ^ — °,AT\ ' (6) 

Jo + mt{—)r 

and parameter a is also expressed as a function of the BL fit parameters in ref. [ 
^ . Note that in the original BL papers [ ^, ^ the slope parameters were evaluated 
at rrit = to, in an approximation that neglected the transverse mass dependence of 
these values. If the fits are done at low pt, such an approximation can be warranted. 
If the temperature gradient effects are expected to be small, the phenonienological 
formula appears in its simplest form, 

T,^To + m{utf. (7) 
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It would be important to experimentally test the transverse mass dependence of 
the slope parameters, i.e. the difference between eqs. (||) and (0). It would also 
be very important to experimentally investigate the rapidity dependence of the 
single particle spectra, and the decrease of the slope parameters in the target and 
projectile fragmentation regions, as suggested by eq. 

Note that the BL hydro also predicts that the rapidity width of the double- 
differential invariant momentum distribution depends on the transverse mass in a 
specific manner, which is very sensitive to the central temperature Tq. In particular, 
the prediction is 

Ay^mt) - Ary^ + ^. (8) 

TTlt 

Hence plotting the squared rapidity width of the single-particle momentum distri- 
bution as a function of would provide a straight line, and its slope parameter 
would yield us an independent handle for the value of the central temperature in the 
hotest and densest regions of the fireballs at RHIC. This information would be very 
valuable when concluding about the temperature gradient effects and establishing 
the significance of a hot, central region with To approximately the critical temper- 
ature of the quark-hadron phase transition, as suggested by the present analysis. 



5.2. Measurements at large rrit 

More precise determination of the model parameters will be possible if new data 
points become available for the effective source sizes (HBT radii) for pions and kaons 
at higher value of the transverse mass of the pair. The important point is that the 
1/ y/ml scaling of the HBT radius parameters is predicted to emerge in the high 
transverse mass limiting case, with small scaling violating terms, that play a larger 
role at small values of the transverse mass. Such a behavior is indeed seen in the BL 
fits in Fig. 1, in particular, the scaling violating terms are rather apparent in the 
low transverse mass Ri data points. If the temperature gradient effects are not so 
important, the center of particle emission for particles with larger transverse masses 
moves more and more out from the axis, hence more and more transverse flow 
effects result in a faster than l/y^roj decrease in the out component. However, if the 
temperature gradient effects are important, all the 3 radius components will follow 
the 1 / Y^TOt scaling law in the large transverse mass domain. A measurement of the 
effective proton source as a function of mt and perhaps proton-deuteron coalescence 
measurements would be also tremendously useful to clarify the status of this scaling 
law. 



5. 3. Non-identical particle correlations 

The BL model could be extended to study non-indentical particle correlations [ ^ , 
that provide very important keys to tell what kind of particles are emitted first from 
the hot and dense decaying fireball. In particular, non-identical particle correlations 
are sensitive to both the temporal and spatial separations of the different kind of 
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emitted particles. Preliminary data from the STAR collaboration indicates, that 
protons are emitted closer to the surface, than kaons, and kaons are emitted closer 
to the surface, than pions. In the BL model, the production of heavier particles is 
more and more focussed to the collision axis, the transverse momentum and mass 
dependence of the center of particle production in the transverse plane is given by 
eq. (132) of ref. [|^, which simplifies at mid-rapidity to 

which implies that at any given value of the transverse momentum pt the center of 
particle production for heavier particles is closer to the collision axis, than that of 
the lighter ones. Note, that the BL-hydro also predicts, that the protons come from 
a smaller effective source than the pions, and due to this effect, some of the pions 
may appear from behind the protons. 

In non-identical particle interferometry, the effects are maximal if the velocity 
of the different particles are approximately the same. It is striking to observe, that 
plotting eq. (|^) as a function of the mean transverse velocity of the particle pair, 
heavier particles with a given velocity are emitted from a more forward region than 
the lighter particles! This feature of the BL model is qualitatively similar to the 
observations by the preliminary STAR data on non-identical particle correlation [ 

This feature of the BL model is closely related to the transverse mass scaling of 
the HBT radius parameters, and it would be very important to check experimentally. 
The BL model predicts a similar, focussing effect in the longitudinal direction too, 
namely that the center of emission of heavier particles at any given value of rapidity 
is closer to the midrapidity, than that of the lighter particles, 

yo-y ,.f.^ 

' Jo 

fz = To sinh(y + 7^), (11) 
tz = TnCOsh{y + r]), (12) 

where y is the rapidity of the particle, ya stands for the value of mid-rapidity and rj is 
the space-time rapidity of the particle in the LCMS {y = 0) frame, while the center 
of particle production is given in the frame of observation. Due to Lorentz boost 
effects, this result implies that within the Buda-Lund picture, heavier particles are 
also emitted earlier, than the lighter ones, if the observation of the temporal and 
spatial sequence is done in the LCMS {y = 0) frame. 



6. Conclusions 

We find that the PHENIX and STAR data on single particle spectra of identified tt~ , 
K~ and p as well as detailed mt dependent HBT radius parameters are consistent 
with the Buda-Lund hydro model as well as with one another. 
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The major difference between the final state of heavy ion collisions at RHIC and 
at CERN SPS seems to be not only the increased freeze-out time and the increased 
transverse flow or Hubble constant and relevantly larger tranverse gemetrical radius 
at RHIC, but, most strikingly, the existence of a hot center located close to the 
beam axis, which evaporates particles with a temperature that is very close to the 
deconfinement temperature, Tq = 202 ± 13 MeV. This effect was not seen at CERN 
SPS, when analysed in terms of the Buda-Lund hydro model. We find that this 
parameter can be determined independently from the temperature and the velocity 
of the surface, for which parameters we found rather conventional values. 

Note also that our findings are in qualitative agreement with numerical results 
found from Humanic's cascade [ |l^, [l^ as well as from URQMD [ It seems 
that at RHIC-1, some of the pions, kaons and protons are emitted directly from a 
rather hot zone, that has (within the errors of the reconstruction) the deconfinement 
temperature of QCD, which even at finite chemical potential is predicted to be below 
the Tc — 172 ± 3 MeV value [ The hot center seems to be surrounded by a 
cooler hadronic matter, with a surface temperature of about Ts w 110 MeV. Clearly, 
precision data from RHIC-2 in a broad rapidity and transverse mass domain are 
needed to finalize the conclusions about the significance of this surprising result, 
in particular, the measurement of the pion and kaon HBT radius parameters at 
TOt ~ 1 GeV would provide a very stringent restriction on the models. 
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Fig. 1. Simultaneous fits to the final PHENIX particle spectra and the final 
PHENIX and STAR HBT radius parameters within the framework of the Buda- 
Lund hydro model. Solid line stands for the best Buda-Lund fits using parameter 
set Fit-I, as given in the first column of Table 1. The quality of the fit is statistically 
acceptable, CL = 29 %. 
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Fig. 2. Simultaneous fits to the final PHENIX particle spectra and the final 
PHENIX and STAR HBT radius parameters within the framework of the Buda- 
Lund hydro model, with Tq = 110 MeV fixed. Solid line stands for the Buda-Lund 
fits using parameter set Fit-Ill, as given in the third column of Table 1. Although 
the fit does not look too bad, statistically the fit is not acceptable, as the confidence 
level of the fit is as low as 10"*^. When comparing with Fig. 1, it is clear that the 
mt dependences of the side and the out radius parameters are not strong enough 
and the number of pions with high transverse momenta is too small. 



18 



T. Csorgo and A. Ster 




Fig. 3. The mean production point of pions and protons with a given velocity 
from the BL hydro model, corresponding to eq. utilizing the best fit values 
given in column I of Tabic 1. Surprizingly, we find that within the BL-hydro model, 
heavier particles are emitted with larger transverse coordinate values, than lighter 
particles, when non-identical particles with the same velocity are compared. The 
BL hydro model predicts according to this figure, that both in the small and the 
large transverse velocity limits the difference between the transverse displacement 
of the pion and proton production points disappears. 



